
ELSEVIER 



Available online at www.sciencedirect.com 



V* ScienceDirect 



Nuclear Physics A 00 (2012) 1-g] 



NUCLEAR 
PHYSICS 



www.elsevier.com/locate/nuclphysa 



Realistic modelling of jets in heavy-ion collisions 

Clint Young 

Department of Physics, McGill University, 3600 University Street, Montreal, Quebec, H3A2T8, Canada 

Bjorn Schenke 

Physics Department, Bldg. 510A, Brookhaven National Laboratory, Upton, NY 11973, USA 

Sangyong Jeon 

Department of Physics, McGill University, 3600 University Street, Montreal, Quebec, H3A 2T8, Canada 

Charles Gale 

Department of Physics, McGill University, 3600 University Street, Montreal, Quebec, H3A 2T8, Canada 



Abstract 

The reconstruction of jets in heavy-ion collisions provides insight into the dynamics of hard partons in media. 
Unlike the spectrum of single hadrons, the spectrum of jets is highly sensitive to q ± , as well as being sensitive to 
partonic energy loss and radiative processes. We use martini, an event generator, to study how finite-temperature 
processes at leading order affect dijets. 



1. Jets in a thermal medium 

Dijets are the aftermath of 2 — > 2 scattering in perturbative QCD: when the scattering is at a momentum scale 
significantly above Aqcd, perturbation theory works well to determine the inclusive cross section. The two outgoing 
partons are born high in virtuality and evolve down through collinear processes; this leads to two narrow cones of 
partons centered on axes determined by the original partons' momenta. Hadronization is non-perturbative and soft, 
and leads to hadrons centered on these axes. Because the initial partons had relatively small total transverse momenta, 
the total transverse momenta in the two jets should be roughly balanced: p rl « -pri- 

This picture leads one to think that the total cross section for dijets can be easily compared with perturbative cal- 
culations for inclusive 2 — > 2 scattering. However, even in proton-proton collisions with relatively small backgrounds, 
a method for deciding what particles belong in a jet must be made. A deceivingly simple method is to use the anti-^7- 
algorithm [1 1: a metric defines a "distance" between two 4-momenta; for infrared and collinear 1 — > 2 processes, 
this metric is defined so that the distance between the two final momenta are small. A list of distances between all 
pairs of momenta in a given event is then determined and ordered from smallest to largest, and the list is shrunken 
by combining these momenta, starting at the beginning of the list and stopping once all distances are above some 
lower limit. The final momenta, when they are sufficiently large, can be called reconstructed jets. Because the shortest 
distances are added first, the algorithm is insensitive to infrared and collinear physics. 

The large cross sections for jet production, as well as the large energy separation between jet energies and the 
highest temperatures reached, allowed jets to be one of the first heavy-ion observables examined at the Large Hadron 
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Figure 1 : (Colour online) The effect of elastic and radiative processes on dijet asymmetry. Left: dijet asymmetry of proton-proton collisions (black) 
at 2.76 TeV compared with the Aj distribution in 0-10% lead-lead collisions considering only elastic processes. Data from ATLAS for p+p and 
0-10% Pb+Pb collisions are shown, with (red and black) solid circles. Right: the dijet asymmetry in 0-10% lead-lead collisions considering only 
radiative processes. 



Collider (2l 0. A large sample of dijets were reconstructed from lead-lead events measured at ATLAS, and the 
distributions in dijet asymmetry Aj = (£> - E < )/(E > + E<), where £> (£<) is the energy of the (second-)largest jet, 
were determined for several centrality classes. In the 0-10% centrality class, the distribution dN/dAj was determined 
to be flat from Aj = up to Aj = 0.5, significantly different from all preceding proton-proton analyses where the 
distribution is strongly peaked at Aj = 0. 

We consider the propagation of the partons of a jet through a deconfining, thermal medium. Determining the 
propagation perturbatively leads first to considering elastic 2 — > 2 scattering at tree-level, and using the hard thermal 
loop (HTL) approximation leads to the rate 

dJ r4 1 = tAto f n7 M W**(P +k-p'- fc>(**)[l ± »(«*)]■ (i) 

When interference between multiple elastic scatters is small, these rates can be convoluted with a distribution of 
partons and calculated in the rest frames of the medium and at temperatures determined by a model for the heavy-ion 
collision. Figure [T] shows how convoluting elastic rates can lead to some dijet asymmetry. Using a s =s 0.3, and only 
elastic processes, does not lead to a large average dijet asymmetry. 

However, elastic scattering is not the only process contributing to the evolution of jets at leading order in g. At 
high energies, 2 — > n scatterings where n > 2 become collinearly enhanced. All that keeps these diagrams from 
being singular in the collinear limit is the thermal mass squared ~ g 2 T 2 , which causes scatters with gluon emission 
to contribute at a smaller order in g than one might naively expect. However, at the high temperatures where HTL 
effective theory applies, the medium is also very dense, meaning that the Landau-Pomeranchuk-Migdal (LPM) effect 
must also be taken into account. Multiple scatterings can be resummed using an integral equation [4|. Figure [T] also 
shows the effect of radiative splittings alone. These radiative splittings, which are in addition to the splittings that 
occur in the evolution from high to low virtuality, lead to a distribution in dijet asymmetry almost the same as the 
distribution predicted for proton-proton collisions at this energy. This should not be a surprise: the anti-fc^ algorithm 
is designed to be insensitive to collinear processes. 

However, what happens when both elastic and radiative processes are taken into account? Radiation evolves the 
distribution of partons in a jet down to small momentum fractions. With elastic scatterings occurring, the physics of 
q also is taken into account. Given a parton of momentum p, an elastic scattering transferring transverse momentum 
~ gT will deflect the parton by 69 ~ gT/p. Thanks to the evolution down to small momentum fractions, many more 
of these partons will have small p and will be deflected more. This identifies a possible source for dijet asymmetry in 
the absence of strong energy loss: enhanced decollimation of the jet caused by in-medium jet evolution. 
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2. The importance of full Monte Carlo 

The evolution of a distribution of partons in a thermal medium can be described with this integral equation: 
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If these rates are small unless k is small compared with p, then one might want to make a "diffusive approximation" 
that turns these equations into differential equations. However, for massless partons, this was demonstrated in Figure 
3 in O to pose difficulties: the initial distribution 6(E - 10 GeV) was evolved for 5 fm/c according both to integral 
equations and the Fokker-Planck equation. While both approaches determined similar mean values, the Fokker-Planck 
equation predicted a significantly smaller distribution at small E than solving the integral equation did. 

Our description of jet modification depends critically on determining the jet shape properly at small momentum 
fractions; it is clear that the full integral equation must be solved. However, solving these equations in a rapidly 
evolving and expanding thermal medium is prohibitive and indeed, undesirable. Monte Carlo simulations are ideal 
for solving these integral equations: they are not significantly complicated by the addition of new processes like other 
numerical solutions, and they are adaptable. 



3. Dijet observables at RHIC and the LHC 

In 0, we showed how dijet asymmetry might be described as the result of parton evolution at finite temperature, 
using martini, an event generator for high-momentum observables in heavy -ion collisions Q. With it, jets were 
sampled using pythia8 (8), and evolved according to finite-temperature rates using the results of music, a solver of 
the 3+1 -dimensional relativistic hydrodynamical equations which describe well the bulk observables at RHIC and the 
LHC J9). After hadronization, jets are reconstructed using fastjet IflOl . to be consistent with the techniques used by 
the experimentalists. 

The main results of our work on dijets are shown in Figure|2] Applying both the elastic and radiative processes to 
the partons in a jet leads to agreement with the results of both ATLAS and CMS, using a s ~ 0.3, unlike the results 
found using elastic or radiative processes alone. 



4. Raa and the effect of running coupling and finite-size effects 

Finally, we present preliminary results for the charged particle nuclear modification factor RpbptApr) in the 0-5% 
centrality class at the LHC, in Figure[3] Three curves are compared with data: first, the results using the same radiative 
rates as was used to describe dijet asymmetry are plotted with the solid (red) curve. A slow rise in Raa with pj is 
observed, consistent with the LPM effect IfTTI . Next, the dashed (green) curve shows results with a parametrization 
of finite-size effects. For a medium with finite size, the interference between vacuum and medium-induced splittings 
suppresses the radiative rates at early times Ifl2l . Finally, the dotted (blue) curve shows the results taking account 
of both finite size effects and the running of the coupling. The momentum scale of the splitting vertex is given by 
(lefj = qtf = qk/ (k%), leading to {\kr\) = (qk)'* . Momentum transfers in elastic collisions are space-like and of order 
gT and their running is ignored. 
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Figure 2: (Colour online) dN/dAj, including both elastic and radiative processes, with experimental cuts, compared with dat from ATLAS and 
CMS. 
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Figure 3: (Colour online) The nuclear modification factor Rpi,pi,(Pt), including the effect of finite-size and running coupling, compared with 
preliminary CMS results. 

At this point, both effects are applied approximately to total radiative rates and not to the differential rates. To 
include the effect of running coupling without this approximation is simply a matter of recalculating the rates. How- 
ever, exact implementation of finite-size effects poses a challenge to Monte Carlo simulations, because they depend 
on the previous positions of the partons and the evolution of the medium. Perhaps using the initial position of a given 
parton, and the temperature of the medium at this point, can be used to estimate the temperature profile for the parton's 
evolution with enough accuracy that the finite-size effect can be taken into account accurately. 
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